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We study the electromagnetic properties of a system that consists of an electron background and 
a neutrino gas that may be moving or at rest, as a whole, relative to the background. The photon 
self-energy for this system is characterized by the usual transverse and longitudinal polarization 
functions, and two additional ones which are the focus of our calculations, that give rise to birefrin- 
gence and anisotropic effects in the photon dispersion relations. Expressions for them are obtained, 
which depend on the neutrino number densities and involve momentum integrals over the electron 
distribution functions, and are valid for any value of the photon momentum and general conditions 
of the electron gas. Those expressions are evaluated explicitly for several special cases and approxi- 
mations which are generally useful in astrophysical and cosmological settings. Besides studying the 
photon dispersion relations, we consider the macroscopic electrodynamic equations for this system, 
which involve the standard dielectric and permeability constants plus two additional ones related to 
the photon self-energy functions. As an illustration, the equations are used to discuss the evolution 
of a magnetic field perturbation in such a medium. This particular phenomena has also been con- 
sidered in a recent work by Semikoz and Sokoloff as a mechanism for the generation of large-scale 
magnetic fields in the Early Universe as a consequence of the neutrino-plasma interactions, and 
allows us to establish contact with a specific application in a well defined context, with a broader 
scope and from a very different point of view. 



I. INTRODUCTION AND SUMMARY 



This work is concerned with the electromagnetic properties of a medium that consists of a matter background, such 
as an electron plasma, and a neutrino gas that moves, as a whole, relative to the matter background. Technically, 
the quantity of interest to us is the photon self-energy, from which the dispersion relations of the photon modes that 
propagate in the medium can be obtained, and from which other macroscopic quantities of physical interest can be 
determined. 

Some aspects of this composite system were studied in Ref. using the methods of real-time finite temperature 
field theory (FTFT)|28l| which, from a modern point of view, provides a natural setting for studying the problems 
related to the propagation of photons in a medium. Largely stimulated by the work of Weldon0,y,yj, a convenient 
technique employed in FTFT is to carry out the calculations in a manifestly covariant form. As is by now familiar, this 
is implemented by introducing the velocity four-vector of the medium, in terms of which the thermal propagators 
are written in a covariant form. In this way, covariance is maintained, but quantities such as the photon self-energy 
depend on the vector it M in addition to the kinematic momentum variables of the problem. Generally, for practical 
purposes the vector is set to (1, 0) in the end, which is equivalent to having carried out the calculation in the rest 
frame of the medium from the beginning, and this is usually the relevant physical situation. 

However, as noted in Ref. [l|, the system we are considering provides an example for a novel application of FTFT. 
A distinctive feature of this system is that the matter background on one hand, and the neutrino gas on the other, 
each is characterized by its own velocity four-vector. Thus, if wc denote by the velocity four-vector of the matter 
background, and by v M the corresponding one for the neutrino gas, then we can take the matter background to be at 
rest, so that 

u" = (i,6), (l.i) 

but we must keep 

v ti = (v°,V). (1.2) 

Therefore, the physical quantities such as the photon self-energy depend on the momentum variables and as usual, 
and in addition on the vector v^. This additional dependence can have physical effects that cannot be produced by 
the stationary background alone. 
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FIG. 1: A class of diagrams that contribute to the photon self-energy in the presence of a neutrino background. The filled 
circle in diagram (a) represents the neutrino electromagnetic vertex that is induced by the interactions with the electron 
background p. Isl Hfl 111 YVA Il3j| . which can be determined by calculating the off-shell amplitude corresponding to diagrams (b) 
and (c)Q. 
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FIG. 2: A class of diagrams that contribute to the photon self-energy in the presence of a neutrino background. The double 
electron line represents the thermal electron propagator, with the effects of the neutrino gas taken into account. 

The focus of attention in Ref. Q was the effect that the collective neutrino-plasma interactions could have on the 
stability of such systems 0, 0, in analogy with the stream instabilities that are familiar in plasma physics research, 
examples of which are discussed in many textbooksQ. The calculation was focused on the contribution to the photon 
self-energy from the set of diagrams represented in Fig. ^ 

On the other hand, it was recognized some time ago |l5l fl6l Il7l fl8| that the presence of neutrinos induces birefrin- 
gence effects in the medium, similar to those that exist in materials that exhibit natural optical activitv^L^. These 
effects are a consequence of a new term, denoted by np in Ref. |l5j| . that shows up in the photon self-energy, which is 
proportional to the neutrino-antincutrino asymmetry and is odd under various discrete space-time symmetries. How- 
ever, the previous calculations of these effects have taken into account only the neutrino background. In particular, for 
the purpose of determining irp in Ref. p"?|. the background was taken to consist only of neutrinos and antineutrinos. 
Therefore, for all the other particles the vacuum propagator was used in the calculation. 

Here we note that, in the presence of neutrinos, there is an important contribution to the photon self-energy which 
is proportional to both the neutrino asymmetry and the electron number density. It is best specified by referring 
to the diagram depicted in Fig. [21 As indicated by the double line in this diagram, the electron propagator to be 
used in the calculation is the dressed propagator with the correct electron self-energy in the medium. The electron 
self-energy includes the effect of the neutrino gas, as represented by the diagrams shown in Fig. |3J As we show, this 
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FIG. 3: One-loop diagrams for the electron self-energy in a neutrino gas. 



contribution to the photon self-energy is of no consequence for the instability issues, but it yields a novel contribution 
to the optical activity of the system, that can lead to interesting and significant physical effects. In particular, we note 
that the contribution to np calculated in Ref. |l7j | is proportional to the photon momentum squared q 2 . In contrast, 
the contribution that is determined here has no such term, and it is the only one that survives in the q — > limit. 
This feature leads to effects that can manifest themselves at a macroscopic level, in the static and long wavelength 
electromagnetic regime. 

In this work we compute the contribution to the photon self-energy that is proportional to both the neutrino- 
antineutrino asymmetry and the electron density, and consider the consequences for the propagation of photons and 
for the electromagnetic properties of the medium. The calculation is based on the application of FTFT to calculate 
the photon self-energy diagram shown in Fig. [21 using the electron propagator that includes the effect of the neutrino- 
electron interactions. The implicit assumption is that, in its own rest frame, the neutrino gas has a momentum 
distribution function that is parametrized in the usual way. The presence of the neutrinos gives rise to two additional 
polarization functions besides the ordinary longitudinal and transverse polarization functions the photon self-energy, 
that we denote by irp and tt' p . A non-zero value of irp by itself leads to optical activity effects, while ir' p induces 
anisotropic effects. Aside from the implications for the photon dispersion relations, the results for the photon self- 
energy can be interpreted in terms of the macroscopic electromagnetic properties of the system. In particular, the 
effects due to np and 7Tp remain finite in the limit q — > 0, and therefore they can manifest themselves in macroscopic 
effects in the long wavelength and static regimes. As a specific application, we consider the evolution of a magnetic 
field perturbation in the system, and we arrive at an equation for the dynamics of the magnetic field that has been 
suggested by Semikoz and Sokoloff 20] , in the context of a mechanism for the generation of large-scale magnetic 
fields in the Early Universe as a consequence of the neutrino-plasma interactions. Thus, besides extending the earlier 
calculations already mentioned, the present work makes contact and complements this recent work, which is based 
on a treatment using the kinetic equations of the neutrino-plasma system. The results and formulas we present are 
applicable in a variety of situations where the neutrino interactions with the other background particles are important, 
and the method we employ could also be useful in the study of similar problems that may arise in other contexts. 

A word about the strate gy o f our calculation is in order. In principle, the diagrams involved are numerous, and many 
of them are shown in Ref. [l7| . However our calculation is simplified for the following reason. Since we are interested 
in the terms that contain both the neutrino and the electron distribution functions, the momentum integrations are 
effectively cutoff. To order l/M^, we can then replace the W, Z boson propagators by their local limit, so that the 
neutrino-electron interactions can be approximated by the local Fermi interactions. In that limit, the set of diagrams 
that are relevant to extract the contribution that we are seeking, to order 1/M^, collapses to the class of diagram 
represented in Fig. [21 which we are considering in the manner we have indicated. 

We begin, in Section[n]by writing down the expression for the electron propagator that will be used in the calculation 
of the photon self-energy, which includes the thermal effects as well as the effects of the neutrino-electron interactions. 
In Section IIIII the self-energy tensor is calculated and in terms of the usual transverse and longitudinal components 
and the two additional components that we have mentioned, and the integral formulas for the latter are obtained and 
evaluated for some cases in Section llVl The photon dispersion relations for this system are considered in Section Ivl 
the application to the study of the macroscopic electromagnetic properties of this system is in Section IVll and Section 
IVIII contains our conclusions. 
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II. ELECTRON PROPAGATOR 



The electron self-energy function S e (p) is determined by computing the diagrams depicted in Fig. [3] A straight- 
forward calculation, using the thermal propagator for the internal neutrino lines, yields to order X/My^ 



(2.1) 



where 



Xv — \/2G F 
\ A = V2G F 



ai ' > \n Ua - n Va ) +n v<: - n v 

a=e,fj,,r 



7 E ( 

( n ^« ~ ~ ( n v e — n vJ 



ct— e,£*,r 



(2.2) 



( z) ( z) 

Here the coefficients a e \b e are the vector and axial vector neutral-current couplings of the electron, 

1 

2 

b (z) = I 
2 ' 



ai z) = -- + 2sm 2 w 



(2.3) 



and the n Va , n Pa denote respectively the total number density of each neutrino or antineutrino specie in the medium. 
The electron propagator to be used in our calculation of the photon self-energy is given by 



where S F — j°S F "f°, and 



with 



S e (p) = S F {p) ~ [S F (p) ~ S F (p)] Ve(p ■ u) , 
Ve(p) = 0(p ■ u)f e (p ■ u) + 6(-p ■ u)f s (-p ■ u) . 



(2.4) 
(2.5) 



fe(x) 



h{x) 



i 



g/3(£C-/i e ) _|_ i 

1 

e /3 e (a;+/j e ) _|_ I 



(2.6) 



Here e is the inverse temperature and /i e the chemical potential of the electron gas. In addition, S F is the electron 
propagator in the presence of the neutrino gas, which is given by 



S F \p) = S^ 1 - E e (p) , 



where Sq is the free propagator in the vacuum 



So 



•p + m e 



p2 _ m 2 _|_ j e 



and E e has been given in Eq. I|2.1|) . Therefore, to leading oder in 1/Myy, 

S F {p) = S (p) + S (p)S e (p)5 (p) • 
Substituting this in Eq. Ij2.4|l . we obtain the electron thermal propagator to be used in our computation 

S e = Sq + St + S' + S'rp , 
where So is given in Eq. (|2.8(l . St is the usual thermal part of the electron propagator 

S T (j>) = 2m8(p 2 - m 2 e )r] e (p • + m e ) , 

while 

(^ + m e )I] e (p)(^ + m e ) 



S'(p) 



(p 2 — m 2 + ie) 2 



(2.7) 

(2.8) 

(2.9) 
(2.10) 
(2.11) 

(2.12) 
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and 

S' T {p) = -2mS'(p 2 - m^e(p • u)(j/> + m e )Ze(p)ty + m e ) . (2.13) 

In Eq. H2.13fl . S' denotes the derivative of the delta function with respect to its argument. 

In some respects, the present calculation resembles the calculation of the photon self-energy in an electron back- 
ground in the presence of a magnetic field |2ll l22j . In that case, the calculation involves the use of the thermal 
generalization of the Schwinger propagator, which takes into account the B field and the thermal effects of the elec- 
tron background. In the present case, the neutrino current acts as the external field, playing the role of the B field in 
the former case. 

III. PHOTON SELF-ENERGY 

We decompose the photon self-energy into various parts according to whether or not they depend on the neutrino 
densities. Therefore, discarding the term that is independent of the particle densities, we write 

7^ = 7$ +7$ +7^, (3.1) 

where is the purely electronic contribution, while ~k^ v > , which depends on the neutrino distributions but not on 
the electron distribution, corresponds to the contribution to the photon self-energy that was computed in Ref. [l7| . 

Although the results for , and of course are known, we state them below in the form that will be useful for later 
reference. On the other hand, ir^" contains the terms that depend on both the neutrino and electron distributions, 
and is the focus of the present paper. To determine it, we start from the expression corresponding to the diagram of 
Fig.0 

/d 4 p 
75^4 Tr h^ s eip + q)luiS e (p)] , (3.2) 

where S e (p) is given by Eq. I|2.4[) . with Sf determined from Eq. I|2.7|) . Substituting Eq. (|2.1UI) into Eq. (|3.2|) and 

singling out the terms that contain the electron distributions as well as the neutrino densities we obtain, 



d 4 p 
(2t7) 

+ Tr llJ iS l {p') lv iS T {p) + Tv ltl iS T {p') lu iS\p)} , (3.3) 



inffi = e 2 / {^l^S [jp')-i v iS' T {p) +Trj ll iS' T (p')j„iS (p) 



where 



p'=p + q. (3.4) 



A. 4 e j 



The pure electron term is given by 



^t )= e 2 / 77^{^l^ s o{p')luiS T {p) J r^l P ,iS T {p')luiS G {p)} . (3.5) 



(27T) 

Computing the traces and carrying out the p° integration, Trffi can be written in the form 



where 



L^u = Zp^pv +p ll q u + q^p v - p ■ qg^ v . (3.7) 
The fact that "nffl is symmetric and satisfies q^irffl — q u ir^u — implies that it is of the form 

tt$ = n^ ) R^ + n^ ) Q fiU , (3.8) 
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where 



9iiv Q(n> ■ 



with 



and 



q^q v 



In general tt^l are functions of the scalar variables 



q ■ u 



Q = v 7 ^ 3 ? 



(3.9) 

(3.10) 
(3.11) 

(3.12) 



which have the interpretation of being the photon energy and momentum, in the rest frame of the electron gas. The 
functions L are determined by projecting Eq. (|3.6J) with the tensors and Q^. This procedure then leads to 



= -2e 2 A 



Q 2 



where 



A e = 



B e 



3o (f £ 



(2tt) 3 2e 

d?p 
(2n) 3 2E 



4e 2 — B 



2m\ ~ 2p ' 1 



(3.13) 



2p-q 



+ (q^ -q) 



(fe + h) 



2(p ■ u) 2 + 2(p ■ u)(q -u) -p - q 



+ (q^ -q) 



(3.14) 



q 2 + 2p ■ q 

While these formulas can be used to evaluate the electronic contribution to the photon self-energy in various situations, 
we need not proceed any further in that direction since, as we have already mentioned, the results are well known 
and we can simply quote the relevant ones when we need them. 



B. 



In addition to the tensors and that are defined in Eq. (|3.9|l . we introduce 



P, 



Q 



Q 



a,.0 



PL = Qe^apq-v 



with 



■< ,f * = v» - v?{u ■ v) . 



Then the results of Ref. can be expressed in the form 

n$=7rP(P^+u-vP tl „), 

where 

M e 2 G F ( q 2 \ 2 

J(q 2 ) is given explicitly in that reference, but its precise value will not be relevant here. We only wish to note the 
presence of the kinematic factor of q 2 in Eq. I|3.18[) . For the purpose of determining the dispersion relations, this 
factor can be set equal to the plasma frequency squar ed uj 2 since in the lowest order q 2 ~ w 2 ,. 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



applications, such the one as we will consider in Scction fVI Bl the appropriate kinematic regime corresponds to taking 
q — > 0, and in those cases 7Tp ' is not relevant. 
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Computing the traces and carrying out the p° integration in Eq. i|3.3[l , 7rj,t^ can be expressed in the form 



where 



with 



and 



fit" = (-4e 2 )(AvT<r> + AjTW) , (3.19) 



nV . / -g^. „) { - Sffl + <1 + (, _ _,)| 

(3.20) 



= ( TO e - P 2 ) [p' M v t> - g^p' ■ v + v,tp' v ] + (2p ■ v)mlg liV + 2p-v [p'^p v + p^ - p ■ p'g^] 

L $ = ( m t-P >2 ) \Pn v v - g^p ■ v + v^p v ] + (2p' ■ v)m 2 e g^ + 2p' ■ v [p'^p v + - p ■ p'g^] 

K ap = (P 2 _ rrfyPaVp + {p 2 + m 2 e )q a Vp - 2p ■ v q a pp , 

K al = ( TO e - P >2 )PaVp + 2m 2 e q a vp - 2p' ■ V q a pp , (3-21) 



d={p + q) 2 -m 2 e . (3.22) 



1. Evaluation ofT^ 

The integral expression for T$ can be simplified as follows. Defining 

Cfj, u = p^v v + v^pv ~ p ■ vg^u , 

D^ u = {mj - p ■ p')g^ v + p^p' y + p'^Pv , (3.23) 

Lffv can be expressed in the form 

= ( m ' ~ [ vXd * D ^ ~ C A + 2 P ■ vD ^ > 

L$ = -dC^ + 2p' -vD^, (3.24) 

where d is given in Eq. I|3.22J) . 9 M = d/dp^ 1 , and we have used the following relation, 

p'^Vv + Vpp'v - p' ■ vg^ u = v x d\D flI/ - C^y , (3.25) 

which can be verified by explicit computation. Using the relations 

{p 2 - ml)6' V - ml) = -5(p 2 ~m 2 e ), 

v*dp5{p 2 -ml) = 2p-vS'(p 2 -m 2 e ), (3.26) 

it then follows that 

L$5'(p 2 - m 2 e ) = v x d x [D^S(p 2 - ml)] - C^S(p 2 - m 2 e ) , (3.27) 

and similarly 

f -D^($-°f, (3.28) 
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where we have used 



1 



d) ~ e <P 
Therefore, from Eq. (j!HU|l . with the help of Eqs. (|3~T7|) and 



d 4 P 
(2tt) 3 



rj e (p ■ u) < -v d\ 



(3.29) 



(3.30) 



which, after a partial integration and then carrying out the integration over p° using the delta function, yields 

3? — «/ ^ {^~ q + <« — "} ■ <«» 

with L M „ given in Eq. I|3.7|l . This has the same structure as the normal electron background contribution, given in 

Eq. I|3.6|l . Therefore, the same arguments used in Section IIII Al can be applied here to conclude that T^u can be 
expressed in the form 



-r 

Q2 I U > ' 



where 



A' e = u-v 
B' = u-v 



d 3 p d(f e + / e ) 
{2n) 3 2E dE 

d 3 P d(f e + f s ) 
{2n) 3 2E dE 



2ml ~ 2p • q 



+ (q^ -q) 



q 2 + 2p ■ q 
2(p ■ u) 2 + 2{p ■ u)(q -u) -p-q 
q 2 + 2p ■ q 



+ (q^ -q) 



(3.32) 



(3.33) 



2. Evaluation of ' T$ 
From Eq. (|3.21|) . using Eq. I|3.26[) we can write 



-K§8'tf-ml) K^5(p 2 -m 2 e ) , 



d 



d 2 



2m e q a V/3 



6(p 2 - ml) S'(p 2 -ml) ] x / q a pp8{p 2 - m 2 ) 
d 2 d +V 9x 



Therefore, 



where 



It 



d 4 p 

(2tt)*) 



ri e {p ■ u) { v x d\ 



p /3 S(p 



2 _ m 2\ 



+ {q^ -q) , 



h = 2mi 



d 4 p 

WW) 



Ve(p-u) 



S(p 2 — ml) S'(p 2 — m 2 ) 
d 2 d 



(9 --9) 



For , by partial integration we obtain 



/; = — v ■ u 



d 3 P g(je + h) p 

{2n) 3 2E 8E P 



1 



q 2 + 2p ■ q 



+ (q^ -q) 



Since the integral is a function only of the vectors q^ and u M , it must be of the form 

Jf = hvP + Itf . 

Substituting it in Eq. (|3~35jl. 

TW=ie„ va0 q a \hv? + hvP\ , 



(3.34) 



(3.35) 



(3.36) 

(3.37) 

(3.38) 
(3.39) 
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and therefore /{ is not relevant. On the other hand, I\ can be determined by using the projection relation 



upll = u h , 



which, from Eq. (|3.37|) . yields the formula 



h 



d 3 P d(f e + h 

(2tt) 3 2£ dE 



q 2 (p -u)-{q- u)(p ■ q) 



q 2 + 2p ■ q 



+ {q^ -«) ■ 



To carry out the integration over p° in Eq. I|3.36[) we write 



(3.40) 
(3.41) 

(3.42) 



and using the usual rule for the integration over the derivative of the delta function after some algebra we obtain 



h = 2mi 



d 3 p d 
{2n) 3 2E dE 



1 ( fe + fe 



2E \q 2 + 2p-q 



+ {q^ -q) 



(3.43) 



D. Summary 

From the results we have obtained, it follows that the photon self-energy can be expressed in the form 

V = TTT-fiV + TTlQ^ + K P P^ U + TTpP^ , (3.44) 

where P^ v an d PL, have been defined in Eq. (|3.15|l . Using Eqs. (|3.8|l and l|3.19[) . together with the results for T$' A ^ 
given in Eqs. H3.32JI and l|3.39|l . the coefficients in Eq. H3.44JI are then given by 



7TT 



< ( A ' e + f- B > e ). 

^+4e 2 \ v (^)B' e , 



Q 



TTp = (u ■ v)lTp' + Tl\ 



where 



(3.45) 



= -4e 2 X A Q[h + (u ■ v)I 2 ] , 



n'^ = ~4e 2 \ A QI 2 , (3.46) 

with 7i j2 given by Eqs. I|3.41|l and (|3.43(l . respectively. 

The terms proportional to Ay give only small correction to the pure electronic terms. For example, consider the 
case of a classical electron distribution, for which we can use 

^f«-/?/e,e". (3.47) 

Then from Eq. (|3~33|) . 

A' e = -pA e , 

B' e = -0B e , (3.48) 

and remembering Eq. I|3.13|l . it follows that the neutrino-dependent contribution to 7Tt,l is smaller than the electron 
term by a factor of /3Ay. Assuming, for illustrative purposes, that the neutrino gas can also be treated classically, 
so that n u ~ T 3 , then (3\y ~ GfT 2 , which is negligible in most situations of interest. Therefore, for all practical 
purposes of interest to us, we will neglect that contribution in Eq. (|3.45|) and use 

0) 

tt l = n[ e) (3.49) 



in what follows. 
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IV. EVALUATION OF tt^T 5 AND tt p 



Useful formulas for 7Tp and 7Tp can be obtained by evaluating the integrals I 12 defined in Eqs. (|3.41|l and 
(|3 .431) in the long wavelength limit. The expressions for 2x,2 in that limit can be obtained by applying the auxiliary 
formula 2d 



d 3 p T{p, q) 
(2tt) 3 [q 2 + 2p ■ q]' 



d 3 p 



-77 Q_ dj 7 nuj *t: 

^ 2 ' dp 2E J ~ 



(2tt) 3 [2Eiu-2p-QY 



and they are valid for 



w,Q< (£), 



(4.1) 



(4.2) 



where {£) denotes a typical average energy of the electrons in the background. In Eq. I|4.1|) the symbol stands for 
the total momentum derivative, 



d d p d 
dp = df> + E~dE' 



Let us consider I2 first. We rewrite Eq. (|3.43|l in the form 



h = 



r 



2uf 



where, to simplify the notation, we have defined 

/ = 
/' = 

Then by direct application of Eq. (|4.1(l we obtain 



d 3 p 

(2tt) 3 1 q 2 + 2q ■ p [q 2 + 2q- p} 2 

fe + fe 



+ (q^ -q) ■ 



(4.3) 



(4.4) 



E 2 ' 

1 d ffe + fe 



EdE 



E 



(4.5) 



h = 



d 3 p 



4 J (2tt) 3 E[luE~p-Q} 2 
d 3 p 1 



p-Q 



df_ 

dE 



2uf 



4 J (2tt) 3 E[ujE -p-Q] 



(4.6) 



Furthermore, denoting by v the average velocity of the particles in the background, the following approximate formulas 
are useful for practical applications, 



m e r d 3 p 1 d 
4 J (27r) a E dE 



h = 



(4.7) 



which are obtained from Eq. I|4.6|l by taking the static limit u> = 0, or the Q = limit, respectively. The remaining 
integration can be carried out simply for particular cases of the distribution functions, and we consider some of them 
below. 

Regarding I\ , it can be written in the form 



2Q- 



d 3 P 



where 



2tt) 3 q 2 + 2p ■ q 

d(fe + fe) 



dE 



-Q 



+ {q^ -q) , 



ujp- Q 
E 



(4.8) 



(4.9) 
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and applying Eq. ^4.1(1 . 



4Q2 7 ( 27r )3 ^ dp- + £ i; ^ g 



To compute the momentum derivative we use 



dF 1 „ 2 d 2 (f e + f s ) _ ~ d 
-M = - Q dE* +0Jp - Q dE 



1 S(f e + fg) 



E dE 



which, using p '■ Q = Elu — (Elu — p ■ Q), we write as 



~dE 



l2 9 2 (/e + / g ) d 
+ U) £j 



dE 2 



dE 



E dE 



lu(Elu — p-Q) 



dE 



1 d(f e + fg) 



E dE 



In addition 



dF 1 = to d(f e + fg) 
dp E dE 

and therefore, remembering Eq. I|4.3[) . 

dFi _ 2 p-Qd 2 (f e + f s ) p-Q 



Q 



E F ^ Q -^W = «- E 



dE 2 



1 0(f e + fg) 



E dE 



When this is substituted in Eq. (|4.1U|) . the second term gives zero by symmetric integration, and we obtain 

t( n\- u ' v q2 f d3p 1 d 2 (/ e + / g ) p-Q 

ll[UJ ^> 4 Q 2 J (2tt) 3 E dE 2 Elu -p-Q' 

To evaluate I± in the Q — ► limit, we first expand the denominator in the integrand, 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



(4.14) 



(4.15) 



Ecu -p-Q Ecu 



P-Q 



Euj 



(4.16) 

The term proportional to a single power of p ■ Q integrates to zero, while in the quadratic term we can put (p ■ Q) 2 — > 



1 p 2 Q 2 . In this way, in analogy with Eq. I|4.7|l . we obtain the approximate formulas 



h 



— J (2^f"e — aw — l w ^ V Q> 

u-v r _rffp p 2 d 2 (U+f e ) , - m 

'— ) J2^JW — aW — (u^vQ) 



(4.17) 



In what follows we give explicit results for several useful cases. 

A. Non-relativistic, non-degenerate gas 

For a classical and non-relativistic gas, using Eq. I|3.47[) and fg ~ 0, Eqs. 1)4. 7|) and l|4.17|l yield 



h = —u ■ v I 2 



with 



&- ^«vQ) 



8 m.; 



(lo > vQ) , 



(4.18) 



(4.19) 



where n e is the total number density of electrons. 
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B. Relativistic, non-degenerate gas 

In this case we also assume for simplicity that f g — f e - Then, for I\, we use Eq. (|3.47() twice, and in the limit 

m e — ► 



u ■ v 



1 (w < vQ) 



h = -jj x \ (4.20) 

For J2 we have to be careful because the integral is not defined for m e = 0. As we show below, the result is 

h = ±, (4.21) 

independently of whether to <C vQ or w 3> vQ. To illustrate how we have proceeded let us consider the u) 3> vQ case 
in some detail. By carrying out the angular integration, followed by an integration by parts, from Eq. (|4.7|l we obtain 
the formula 

By making the change of variable p = m e tan 8, 

J d P^fi = ^2 / ^(/ e + /e)cOS0 

= 4t [ d9e-P maece cos9, (4.23) 

which we evaluate by using the Taylor expansion of the exponential. In this way we obtain 

f 4- f- 2 

dP^TT 1 -— [l + OKO], (4-24) 

E 3 TOg 

and substituting it in Eq. 14.22|) we arrive at Eq. I|4.21f) . In similar fashion, for uj <C vQ, 



8tt 2 J ^ EdE \ E 
Using Eq. (|3.47|) and then making the change of variables as above, we obtain 



which leads to Eq. (|4.21() . 



♦sb^-sI 1 ^' (4 - 26) 



C. Degenerate gas 



For a degenerate gas, whether it is relativistic or not, 

f e = 6{E F -E), 

h - 0, (4.27) 



where Ep — \/pp + m e 1S the Fermi energy, with 

PF = (3vr 2 n e ) 1/3 (4.28) 
being the Fermi momentum. Then we obtain in this case 



E 



1 /u-v\ 



(w « vQ) 



E F 



1 _ 2 Pf 

3Ej 



(lu > vQ) , 



^ (w<i;Q) 



1 \ I PF 



h = ( ^ ) X { (4.29) 
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Furthermore, in the non-relativistic (NR) or the extremely relativistic (ER) case, this yields 

ANR) _ T (NR) _ (U-V\ \ VF 



8tt 2 



ir~> = -u-vl?">=- X (4.30) 



AER) _ 
1 1 ~ 



and 



independently of whether u> <C orw> uQ. 



. 1 (w« vQ) 

| (w > vQ) , 



if K) = JL (4.32) 



D. Explicit formulas for 7Tp"' ) and 7rp e "' 

The formulas that we have obtained for I\ i allow us to evaluate 7Tp V ^ and Trp e ^ for typical situations of interest 
in physical applications. We specifically consider the regime 

w>wQ, (4.33) 

which will be of special interest to us in Section IYl Bl Thus, for example, if the electron gas is degenerate, 



O) _ / e 2 \ X A p 



9 3 

r l(eu) _ _ ( X A PF 

2ttV E f 



= -l^)^&Q- (4-34) 



These expressions hold whether the electrons are relativistic or not. On the other hand, for a classical and relativistic 
electron gas, 

(eu) ( 2e 2 



TT p = - I — — ) X A Qu-V 



3ir 2 

2 ' 

' XaQ ■ (4.35) 



Regarding the value of X A , from Eq. (|2.2[l . 

Aa = —= [{n Vr ~ np T ) + (n Uft - n p J - (n„ e - n p J] . (4.36) 
Denoting by \x Vu . the chemical potential of the neutrino of flavor x = e, fi, r, the approximate formulas 

X A = ^=T 3 (& T + - U ) (M,, « T) (4.37) 

or 

= (m" t + < ~ /£. ) (P*. » T) , (4.38) 

can be useful. 
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E. Discussion 



It should be noted that the contributions TVp and ir'p that we have calculated in this work, have a very different 

kinematic dependence on the photon momentum if we compare them with the term itp ' that was calculated in 

Ref. 01 an d the analogous quantity calculated in Ref. ^^|. In particular, Kp\ which does not depend explicitly on 
the electron distribution, is proportional to q 2 , as we have indicated in Eq. I|3.18[l . For the purpose of determining 

the photon dispersion relations, since q 2 is of the order of the plasma frequency squared, the value of 7Tp ' turns out 

to be proportional to the electron density and in fact comparable to the values of TTp^ and ir'p^ . In the absence of 

the electrons, ■Kp' and ir'p^ are of course zero and, since q 2 ~ 0, -Kp^ is negligible. Thus, in a pure neutrino gas 
the term calculated in Ref. [li|. although it is of higher order in l/M^, is the dominant one. On the other hand, 
in applications such as the one that we consider in Section IYl Bl in which the relevant kinematic limit is q — > 0, 
corresponding to the long wavelength and static regime, Wp^ does not contribute and iTp"^ and ~Kp V ^ are the only 
relevant ones. 



V. DISPERSION RELATIONS 

In the presence of an external current j'ySi , the electromagnetic potential in the medium is determined from the 
field equation, which in momentum space is 

[-q% v +^ v ]A v =j^ xt) . (5.1) 

The photon dispersion relations are determined by finding the solutions of the homogeneous equation, i.e., with 
j(P xt ) _ o Xo specify the various modes it is convenient to introduce the basis vectors e^ 2 3 , which we defined as 
follows. 

For a given photon momentum Q, we define the unit vectors (i = 1, 2, 3) by writing 

Q = Qe 3 , (5.2) 

with ei t 2 chosen such that 

ei • S3 = ei • e 3 = , 

e 2 = e 3 x ei . (5.3) 

For the problem that we are considering, without loss of generality, we can choose the vectors ei ( 2 such that V lies in 
the 1,3 plane. Thus, the unit vector 



has the decomposition 



where 



The vectors 6j 23 are then defined by 



V = £ , (5.4) 
V = cos6*e 3 + sin#ei . (5.5) 
cos6 = Q-V. (5.6) 



c? - (0,ei) ; 
4 = (0,e 2 ). 



1 



3 - -_(Q >w g a ), (5.7) 
which form a basis of vectors orthogonal to <? M . They satisfy the relations 
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from which we can immediately read off the matrix elements of the tensors R,Q, P between any pair of the basis 
vectors 2 3 . In addition, the only non-zero matrix elements of the tensor P' are given by 



f^f v P' 



Q 



Vsind, 



Ji v pi _ _ e P e v p' 
t l t 2 Jr nu ~ t 2 t l Jr nu 



Q 



Vcost 



To find the propagating modes, we express the polarization vectors in the form 
where the coefficients at and the corresponding dispersion relations are to be found by solving the equation 

[-q 2 g^ + £' y = • 

With the help of the relations given in Eqs. (|5.8|l and l|5.9|l . this equation can be written in matrix notation 

(q 2 -n)a = 0, 



where 



a i 

«3 



and 



/ 7TT 

II = «7Tp — i^Vir'p cos ( 

V 



-i-Kp + i^Vn'p cos 9 



Tlx 



Q 



i^-Vn'p sin i 



7TL 



(5.9) 

(5.10) 

(5.11) 
(5.12) 

(5.13) 



(5.14) 



In what follows we consider some particular cases of this equation, whose solutions reveal the structure and the main 
features of more general ones. 



A. V = case 

In this case the n' P term in Eq. (|3.44|) is absent, and consequently the dispersion relations have the same form as 
those studied earlier in Ref . [l5] . The longitudinal mode has the dispersion relation 



lo 2 - Q 2 - n L = , 



and corresponding polarization vector 



e L = 4- 

The polarization vectors of the transverse modes are then given by 

with the corresponding dispersion relations given by 

uj 2 - Q 2 - (tt t ± 7T P ) = , 

respectively. 



(5.15) 



(5.16) 



(5.17) 



(5.18) 
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B. V i= case 

The presence of the ir' P term in Eq. (|3.44|) in general modifies the picture described above, and can give rise to 
new effects. This is due to the fact that the existence of the velocity vector V of the neutrino gas breaks the three- 
dimensional isotropy of the medium, for example not too differently from the way in which it would be broken by the 
presence of an external magnetic field. We consider two special situations. 



1. Propagation parallel to V 

In this case s'm9 = 0, and 

7T T -in P + i^Vir'p 

II = ( in P - i%Vir p n T | . (5.19) 

n L 

Therefore, the longitudinal mode is unaffected, while for the transverse modes the polarization vectors are the same 
as in Eq. (|5.17|) but the dispersion relations are now given by 



- ^7T T + 7Tp - — Vn'pj = 0, 
uj 2 - Q 2 - (n T ~ vr P + ^rVn'p ) =0, (5.20) 



for £±, respectively. 



2. Propagation perpendicular to V 

In this case, setting 6 = ir/2, 







— ITTp 


o \ 


n = 


inp 


Ti'T 


i^fVTT'p 




K 


Q P 


7TL / 



Therefore the modes are neither purely longitudinal nor transverse to Q. Finding the general solution in this case is 
a cumbersome process. In some circumstances it may be appropriate to seek approximate formulas for the dispersion 
relations and polarization vectors as a perturbative expansion in Vtt'p, and that can be extended to other values of 9 
as well. 

Some applications of the optical activity effects induced by neutrinos were considered in Ref. [13. As a rule, the 
effects tend to be small. Our main intention in this section was to indicate how the photon dispersion relations 
are modified by the anisotropic effects produced by a non-zero velocity of the neutrino gas relative to the electron 
background. It can be of interest to carry this further to study the implications of these effects in the context of 
the specific applications considered in Ref. ^7]]> or similar ones. However, that is outside the scope and focus of the 
applications that we have considered, to which we now turn our attention. 



VI. MACROSCOPIC ELECTRODYNAMICS 



Besides modifying the dispersion relations of the propagating modes, the presence of the neutrino gas influence the 
electromagnetic properties of the system in the static and long wavelength regime. To study them, it is useful to 
formulate the results or our calculations using the language of macroscopic electrodynamics. 

In what follows we will assume that the neutrino gas is at rest with respect to the electron gas, that is 



v 1 " = , 



(6.1) 
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since this case already brings out the essential consequences of the presence of the neutrino gas. In this case, P' = 0, 
and the photon self-energy takes the form 

= 7r£ * Rftu + ^hQliv + K P Py, v . (6.2) 

As has been discussed previously |l5l fl6j]. this is indeed the most general of the photon self-energy in an isotropic 
medium, which is the case if V = 0. 



A. Dielectric function 

Introducing the electromagnetic field 

Fnu = -i(qnA v - A^qu) , (6.3) 
the equation of motion, Eq. I|5.1fl . can be written in the form 

-iq^ u =j^+j< ind \ (6.4) 

where 

J [ ; nd) = -tt^A" . (6.5) 

In fact, Eq. (|6.4|) is equivalent to the Maxwell equations, with jj? interpreted as the induced current. For example, 
take the component of Eq. I|6.4fl corresponding to the index v being a spatial index. With the usual definition of the 
fields, 

E = iuoA — iQA° , 

B = iQ x A, (6.6) 



which are related by 



B=-QxE, (6.7) 



the equation is just 

iQxB + iujE = + j( md ) . (6.8) 

Moreover, using Eq. (|6.6J) in Eq. (|6.5J) . the induced current is given in terms of the fields by 

(1 - ei)St + (1 - e t )E t + ie p Q x M\ , (6.9) 



where 







et = 


^? 










£; = 






TTp 


e P = 


LO 2 



(6.10) 



and the longitudinal and transverse components of the electric field are defined by 



Ei - Q(Q ■ E) , 

Et = E-Ei. (6.11) 

Eq. H6.9J) is the most general form of the induced current, involving terms that are linear in the field, and subject 
only to the assumption of isotropy. The quantities e^/ in Eq. I|6.10|l are transverse and longitudinal components of the 
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dielectric constant of the medium. Alternatively, instead of et,i, the dielectric and magnetic permeability functions 
e, \x are introduced by writing the induced current in the equivalent form 



j(ind) _ j 



I -< )E [l-~jQxB + i^tpB 



(6.12) 



where we have used Eq. (|6.7|l . Comparing Eqs. (|6.9|l and (|6.12() . the relations 



1 



(6.13) 



then follow. 

These equations can of course be used to discuss the dispersion relations of the propagating modes and related 
effects. However we do not proceed any further in this direction since that would essentially reproduce what we have 
already considered in Section Instead we turn our attention to another kind of effect that can arise due to the 
presence of the neutrino gas, which can be described on the basis of these equations together with the results we have 
obtained. 



B. Evolution of magnetic fields 

Here we consider the evolution of an initial magnetic perturbation in this medium. In the absence of any external 
sources, the equation for the magnetic field, Eq. I]6.8fl. is 

iQ x B + iwE = ~] {md) , (6.14) 
which using Eq. (|6.12|) we can write in the form 

eujE+ -Q x B + i-fB = 0, (6.15) 
A 4 

with 

7 = -q* = «q. (6-16) 

Since we are interested in following the evolution of B, we eliminate E from this equation by taking the cross product 
with Q and then using Eq. H6.7(l . which yields 

eui 2 B - —Q 2 B + ijQ x B = 0. (6.17) 

We obtain the corresponding equation in coordinate space by taking the long wavelength limit, to 3> vQ, and making 
the quasi-static approximation, lo — > 0. In this limit 

ia 

e^l + — , (6.18) 
where a is the conductivity of the medium. Therefore the equation becomes 

iuxrB - -Q 2 B + ijQ xB = 0, (6.19) 
A* 

or, in coordinate space, 

(j^ = -V 2 B+7VxB. (6.20) 
ot [l 



Although we are using the same symbols \i and 7, in Eq. I|6.20|l they stand for the corresponding quantities evaluated 

o; 

p 



in the long wavelength and static limit, as indicated above. In this limit, np ' does not contribute to 7 in Eq. I|6.16|l . 
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as can be seen from Eq. I|3.18[l . Then using the results for 7Tp ^ summarized in Sectior lTV Dl we can readily determine 



7 for some particular cases. For example, for the relativistic and non degenerate electron gas 



2e 2 X A 

7 = 

while for a degenerate electron gas, 



3tt 2 



e 2 X A ( Pf 



3tt 2 \E f ^ 

with Xa given in Eq. (|4.36l) . For more general cases, 7 can be computed from the formula 



7 = -e 2 X 



d 3 P p 2 d 2 (f e + f s ) , 2 r d 3 P 1 d ff e + f, 



rn 



(2tt) 3 3£ 3 dE 2 e J (2tt) 3 EdE\ E 3 



(6.21) 



(6.22) 



(6.23) 



which follows from the expression for tt p given in Eq. (|3.46|) . together with the formulas for I12 given in Eqs. 14.71) 
and for lo > vQ. 

An equation of the same form as Eq. I|6.20|l was obtained by Semikoz and Sokoloff 2fjj by very different means, 
in their work suggesting a new mechanism for the generation of large-scale magnetic fields in the Early Universe 
as a consequence of the neutrino-plasma interactions. As emphasized in that reference, the equation describes the 
self-excitation of an almost constant magnetic perturbation. 

Our approach is useful in two complementary ways. On one hand, it sheds light on the physical origin of this 
mechanism. Within our formulation, it is a consequence of the optical activity induced by the interaction of neutrinos 
with the other background particles. On the other hand, it puts this mechanism on a firm footing from a computational 
point of view. The optical activity induced by the neutrinos can be characterized by the presence of the irp (and tt' p ) 
terms in the photon self-energy, which also has a definite interpretation in terms of the components of the dielectric 
function that enter in the macroscopic (Maxwell) equations in the medium. Therefore, by focusing on these quantities, 
we are able to give well defined formulas for the parameters that are relevant to this effect, in a way that are applicable 
to a variety of astrophysical and cosmological situations in which the presence of neutrinos is influential. Our work also 
paves the way to incorporate some corrections that can be important in specific applications, such as the anisotropic 
effects that can arise if the neutrino gas has a non-zero velocity relative to the electron background. Further studies 
along these lines are the subject of current work. 



VII. CONCLUSIONS 



We have studied the electromagnetic properties of an electron background, that contains a neutrino gas which is 
either at rest or moving as a whole relative to the background. Apart from the well known longitudinal and transverse 
polarization functions of the photon in a medium, the presence of the neutrinos gives rise to two additional polarization 
functions, that we denote by ttp and tt' p . We have computed that particular contribution to these two functions, n p 

and ir'p^ , that depends on the neutrino-antineutrino asymmetry in the medium as well as the momentum integral 
of the electron (and positron) distribution function. The integrals were evaluated for various specific conditions of 
the electron gas, and explicit formulas that are useful in many situations were given. One of the consequences of a 
non-zero value of 7rp and tt' p is to give rise to birefringence and anisotropic effects in the propagation of a photon 
through that medium. We analyzed various particular situations to indicate how the anisotropics due to the non- 
zero velocity of the neutrino gas can affect the optical activity of the system. A non-zero value of np and ir' P also 
has consequences related to the electromagnetic properties of the system at a macroscopic level, and we considered 
specifically the evolution of a macroscopic magnetic field in this system. We arrived at an equation for the dynamics 
of the magnetic field, that had been suggested in Ref. [2(j as a mechanism for the generation of large-scale magnetic 
fields in the Early Universe as a consequence of the neutrino-plasma interactions. In this way we established contact 
between our work and this particular kind of application, which has been of recent interest. The approach we followed, 
which has been based on the application of Finite Temperature Field Theory, as well as the calculations and results 
that are presented here, helps to put this subject on a firm footing and to set a basis for carrying out further studies 
and applications along these lines using powerful calculational techniques. 
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